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Consider a general dynamical system in n dimensions as follows 

x = f(x,fi), x G M n with n > 1, (1) 

where \i G W is a vector of parameters. Our main objective is to find conditions on 
the parameters such that the system above possesses a first integral. In this paper, we 
concentrate on finding integrals of Lotka-Volterra systems with constant terms. 

The Lotka-Volterra system (LV) has been the subject of intensive studies during 
the past century. The interaction of two species in an ecosystem [1] , a metamorphosis 
of turbulence in plasma physics [2], hydrodynamic equations [3], autocatalytic chemical 
reactions [I] and many more, are of Lotka-Volterra type. Nevertheless, the dynamics 
of such systems is far from being understood. Finding first integrals of LV systems, or 
any dynamical system, gives global information about the long-term behaviour of such 
systems. 

In two-dimensional systems, the existence of a first integral implies that the system 
is completely integrable because the phase portraits are completely characterized. For 
three-dimensional systems, the existence of a first integral means that there cannot be 
chaotic motion as the solutions will live inside the level sets of such an integral function. 
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Here we include constant terms to generalize LV systems. The constant term can be 
considered as a constant rate harvesting. Dynamics and bifurcation analysis of such a 
system have been studied in [5]. 

Many different methods have been developed to study the existence of first integrals 
of LV systems. Perhaps one of the earliest attempts to study the existence of first 
integrals was published by Cairo et al. [6J, who studied the integrability of n-dimensional 
Lotka-Volterra equations using the Carleman embedding method. They sought an 
invariant that may be time-dependent. Cairo and Llibre [7] used a polynomial inverse 
integrating factor to find a condition for the existence of the first integral. The 
Darboux method that uses the relationship between algebraic curves and integrability 
of differential equations has been introduced by Cairo and Llibre |8] to study two- 
dimensional LV systems. Cairo et al. [9] used the same method to search for a first 
integral of two-dimensional quadratic systems. 

The Darboux method has also been used to derive an integral for three-dimensional 
LV systems [lOJ and for the so-called ABC systems, which correspond to particular cases 
of three-dimensional LV. The ABC systems were among the first three-dimensional 
models that were investigated. Grammaticos et al. [H] derived first integrals 
using the Frobenius Integrabilty Theorem method (first introduced by Strelcyn and 
Wojciechowski p2]). Ollagnier [13] has found polynomial first integrals of the ABC 
system. Gao and Liu [H] presented a method that basically relies on changing variables 
to transform three-dimensional LV systems to two-dimensional ones. The existence 
of first integrals follows from integrating the two-dimensional systems. Gao [15] used 
a direct integration method to find first integrals of three-dimensional Lotka-Volterra 
systems. A new algorithm presented by Gonzalez-Gascon and Peralta Salas [16] also 
used three-dimensional LV systems as a test case. 

The approach closest to the work presented here uses the idea of associating a 
Hamiltonian to a first integral of a vector field. It was introduced by Nutku [17J. A 
generalization of this idea to two-dimensional vector fields having a first integral was 
provided by Cairo and Feix [18]. They showed that, through time rescaling, the first 
integral can be considered as a Hamiltonian. Subsequently, Cairo et al. [19] and Hua 
et al. [20] used an Ansatz for their Hamiltonian functions. They assumed that a first 
integral (or an invariant) if is a product of functions, H = P(x, y)(Q(x, y))^(R(x, y)) u , 
where P,Q,R are first degree polynomials, and derived conditions for two-dimensional 
quadratic systems to have a first integral. 

Another Hamiltonian method that has been used works as follows. A general 
system (TjQ) is said to have a Hamiltonian structure if and only if it can be written 
as x = f{x) = S(x)VH(x), where S is a skew-symmetric matrix and if is a smooth 
function. The matrix function S must satisfy the Jacobi identity [21]. Plank [22] has 
used this property to find a Hamiltonian function for two-dimensional LV systems, while 
Gao [21] , using the same property, has derived conditions for three-dimensional systems 
not to be chaotic. 
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In this paper, we will not impose the Jacobi identity on the matrix S since we only 
want H to be a first integral. Thus, it is sufficient to ensure that f{x) can be written 
as S(x)VH(x). Consider the following proposition. 

Proposition 1.1 (McLachlan et al. 1999 [23]). Let f G C r (M n ,M n ) ; r > 1, n > 1, be a 

vector field and H e C (R n , M) be a first integral of the vector field f (i.e. f ■ S7H = 0) for 
allx. Then there is a skew- symmetric matrix function S(x) on the domain {x : VH ^ 0} 
such that f = SVH . 

As a consequence, there is also a skew symmetric matrix function T(x) on the 
domain {x : / ^ 0} such that VH = Tf. We are going to use this idea to find 
first integrals and the associated constraints on the parameters for two- and three- 
dimensional LV systems with constant terms. We call the matrix T an integrating 
factor matrix if and only if 



Making an Ansatz concerning the integrating factor matrix T, we obtain both an integral 
and the conditions on the parameters for its existence. The meaning of curl depends on 
the dimension of the system. For two dimensional systems, curl(T/) is the scalar 



where (TF)i is the i — th component of the vector Tf. In three-dimensional systems, 
curl(T/) = V x Tf as usual. 

In Sees. [2] and [3] we will discuss the results of applying the integrating factor 
matrix approach to two- and three-dimensional Lotka-Volterra systems with constant 
terms, respectively. In either case , we will impose some conditions on the integrating 
factor matrix. We will assume, that its entries can be written as a product of functions 
of a single variable. In the three-dimensional case, we use two different forms of the 
matrix. These assumptions are restrictive, and we do not claim to present the most 
general class of first integrals which can be found with the proposed method. However, 
the most general form of the integrating factor matrix will lead to conditions in the form 
of coupled, nonlinear partial differential equations. The analysis of these equations is 
work in progress. We emphasize, that with the simplifying conditions used in the current 
paper, we reproduce many known first integrals and uncover several new ones both in 
two and three dimensions Each section is concluded with a comparison to earlier results 
from the literature. 

2. Two-dimensional LV systems with constant terms 

In this section, we consider integrals of the two-dimensional Lotka-Volterra system with 
constant terms: 



curl(T/) = 0. 



(2) 




0. 



x\ = fi(xi,x 2 ) = xi(bi + a n xi + a 12 x 2 ) + e±, 
x 2 = h{xu x 2) = x 2 (b 2 + a 21 xi + a 22 x 2 ) + e 2 , 



(3) 



An integrating factor matrix method to find first integrals 4 

where b iy = 1, 2) are arbitrary parameters and ei, e 2 are the constant terms. We 

choose an integrating factor matrix as follows, 

nx 1 ,x 2 )=(^° R -* V (4) 

where the function R = R(xi,X2) is to be determined later on. The matrix T(xi,x 2 ) 
is an integrating factor if and only if the curl of Tf is zero, where / = (fi,f 2 )- As 
is mentioned before, in the two dimensional case, this condition is equivalent to the 
following, 

d{Rh) , d(Rf 2 ) 

The associated first integral H is given by 

H(x 1 ,x 2 ) = J R(x 1 ,x 2 )fi(x 1 ,x 2 ) dx 2 + h{xi), (6) 

where h(x\) is found by imposing dH/dx\ = —Rf 2 . 

Let us assume that R is separable, i.e. R(xi,x 2 ) = A(xi)B(x 2 ). We have the 
following lemma. 

Lemma 2.1. Let fi and f 2 be the functions given in the vector field (Gj) ; then the 
condition (TJJ) determines the forms of A(xi) and B(x 2 ) to be, 



— )x/ /ai \ B(x 2 ) = exp(-^W„T/a 2 
' a 12 a 

where a, (3 and 7 are constants to be determined later on. 



A( Xl ) = exp(=^)x 1 /3 / ai2 , B{x 2 ) = exp(— x -^)x 2 ^ , (7) 
012 a 2 i 



Proof. We substitute R into ([5]) to obtain 

1 dA 1 dB . . . 

012^1X2 + ^3—021^1X2 + F{xi) + &{x 2 ) = 0, (8) 
A dxi B dx 2 

where 

F(xi) = ^-(61X1 + anxi 2 + ei) + 61 + 2ana;i + a 2 \X\, 

G(x 2 ) = —(b 2 x 2 + a 22 x 2 + e 2 ) + &2 + 2a 2 2X 2 + 012^2- (9) 
B 

Applying d 2 j[dx\dx 2 ) to both sides of equation flH]), we obtain a separable differential 
equation in terms of A(x\) and B(x 2 ), which can be solved explicitly. □ 

We then substitute A{x\) and B(x 2 ) to the equation ([8]) to obtain 

7x1 + + /3x 2 + G(x 2 ) = 0, (10) 

or, 

7 x 1 + F(x 1 ) = C, /3x 2 + G(x 2 ) = -C, (11) 
where £ is a constant. Finally, we obtain the following set of equations: 

( 1 h&i) + (7H 1 h 2an + a 2 i)aci H £1 H = C, (12) 

012 012 012 0,12 Q-12 Q12X1 
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and respectively, 



«e 2 , 762 , t_ s , fo ceb 2 , 7a 2 2 , n N __ aa 22 2 , 7 e 2 



+ — + 62) + W ~ + — + 2a 22 + a 12 )x 2 -x 2 2 + = -(.(13) 

a 2 i a 2 i a 2 i a 2 i a 2 i a- 2 ix 2 

The expression above is satisfied if and only if the parameters satisfy the following 

conditions: 

q _ _ &«22 

O12 °21 



/3e x 7e 2 



(15) 



&12 a 21 

ah /3a n 

= 7 H 1 h 2a n + a 2 i (16) 

ai2 012 

a& 2 7a 22 , _v 

= /? 1 h 2a 22 + a 12 (17) 

a 21 a 2 \ 

aei ae 2 763 

1 h 61 = 6 2 . (18) 

Ol 2 Ol2 a 21 °21 

Let us introduce l\ := {3/a%2 + 1 and Z 2 := 7/021 + 1- When a = 0, the system (I15H18I) 
is equivalent to the following system of equation: 

e 1 l 1 = e 1 , (19) 
e 2 / 2 = e 2 , (20) 
han + / 2 a 2 i = — an, (21) 

Ml2 + ^0-22 = -«22, (22) 

hh + 6 2 / 2 = 0, (23) 

which is more comfortable to work on and closely related to Plank [22]. As a 
consequence, in the case of a = 0, our problem can be divided into three different 
subcases, depending on the value of the constant terms. In the following we give the 
corresponding conditions along with the resulting integrals for the cases ei,e 2 7^ 0, 
ei 7^ 0, e 2 = and ei = e 2 = separately. The case e\ = 0, e 2 7^ follows by symmetry 
considerations. We note that the case corresponding to the original Lotka-Volterra 
system where e\ = e 2 = 0, has been discussed by various people. However we have 
included a discussion of this case including some previously unknown special solutions. 
The case where a 7^ is discussed in the subsection (2.4). 

To start with, let us consider ( 11911231) as an overdetermined linear system: 

Al = r, (24) 

where I = (h,h) arid the matrix A and the vector r are to be determined later. The 
system has a solution only if the vector r is orthogonal to the left null space of the 
matrix A. 

2.1. The case ei,e 2 7^ 

We have the case where both the constant terms ex and e 2 are non-zero. Then by (fT9l) 
and ( f20l) . this implies that h = l 2 = 1. Moreover, we can simplify the other conditions 
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to: 

h + b 2 = 0, 2a n + a 2 i = 0, and a 12 + 2a 22 = 0. (25) 

If the Lotka-Volterra system ([31) with non-zero constant terms e and / satisfies conditions 
(|25|) . then it has a first integral that is given by: 

H = bi%ix 2 + a n x 1 2 x 2 - a 22 xix 2 + e x x 2 - e 2 x x - (26) 

2.2. The case e 1 ^ 0, e 2 = 

One of the constant terms, e\ is not zero. This means that the free parameter l\ must be 
1 by (|T9j) . We now have a linear system like (12^1) with a 3 by 1 matrix A and a vector r 
in IR 3 with only one unknown l 2 . Without loss of generality, we assume that the matrix 
A is of rank 1, when A has rank zero we have a trivial integral H = x 2 since x 2 = 0. 
Using the fact that this system must be solvable, we can again find the conditions for 
the existence of the first integral. As we assume that a 21 7^ 0, the solvability conditions 
are given by: 

20-11(322 - a-22 - ai2 = 0, and 2b 2 a u - b x a 2X = 0, (27) 
and we have l 2 = — 2an/a 21 . If l 2 is not zero then the first integral is given by: 

H = x 2 2 (-b 2 xi - ^-Xi - a 22 x\_x 2 + ^). (28) 

2 <2 

However, in the case where the exponent l 2 is zero, the integral is given by: 

H = -b 2 xi - ^-x x 2 - a 22 x\x 2 + ei In \x 2 \. (29) 

Note, that expression fl28l) is also obtained as a limit of result (126]) but the conditions 
stated in the case e± 7^ and e 2 = are more general than those obtained in the limit. 

2.3. The case e\ = and e 2 = 

Finally, if e\ = e 2 = equations (fT9|) and (|20|) are trivial and we have a linear system 
of the form (|2^|) with a 3 by 2 matrix A and a vector e in M. 3 from (I21H23I) . 

If A is of maximal rank then the solvability condition of the linear system fl2^|) is 
given by : 

bia 22 (a 21 - a u ) + b 2 a u (a 12 - a 22 ) = 0. (30) 

If our system satisfies this condition then 

h= a 22 (a 21 -an)_ ^ fa= (31) 

a ll a 22 — G 12 a 21 a ll°22 — a 12 a 21 

When neither Z x nor Z 2 are zero, the first integral is given by: 

i/ = x 1 V 2 ^ + ^i-^ 2 ). (32) 

*2 f2 'l 
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However, either li or l 2 may be zero and if l\ — 0, l 2 ^ and l 2 7^ — 1 then we have 
b 2 — a 22 — 0, On 7^ a 2 i and an integral that is given by: 

rj l 2 M Oil «12 \ /qq\ 

H = x 2 2 {— + — Xl + ,^ x 2 ). (33) 

But when Zi = and Z2 — — 1> this implies b 2 = and 021 = an. It follows that the 
integral is given by: 

H = a 12 In \x 2 \ — a 22 In IseiI — a n — . (34) 

x 2 x 2 

Finally, when both li and l 2 are zero, which implies an = a 22 = 0, the first integral is 
given by: 

H = bi In \x 2 \ + 012X2 — & 2 ln |xi| — a 2 \Xi. (35) 

We remark that the case when l 2 = 0, Zi 7^ 0, /1 7^ —1 and the case when l 2 — 0, l\ — —1 
follow by symmetry considerations. 

Now let us assume that A has rank one. Excluding the trivial case in which one 
column vector is zero, we assume that 

(an, 012, h) T = A(a 2 i, a 22 , b 2 ) T , (36) 

for some A 7^ 0. This leads to the following integral: 

H = xix 2 x , where A = an/a 2 i. (37) 

2.4- The case where a 7^ 

This implies a n = a 2 2 = 0, due to conditions (JHJ). Consider conditions (JTB]) and (fT7|) . 
Eliminating a in both equations gives: 

((5 + ai 2 )6ia 2 i + (7 + a 2 i)b 2 a 12 = 0. (38) 
Consider now condition ( {TBI . After some algebraic manipulation we get: 

(/3 + ai 2 )bia 2 i + (7 + a 2 i)b 2 ai 2 = a(e 2 a 12 - eia 2 i), (39) 
and due to the fact that a 7^ 0, we have 

e 2 ai2 - eia 2 i = 0. (40) 
Finally conditions ( fT5j) gives the following three sub-cases: 

(a) 7 = f3 = 0, (b) 7 7^ 0, /3 = 0, and (c) 7, /3 7^ (41) 

The case where /3 7^ and 7 = follows by symmetry. Note that the conditions in 
subcases (b) and (c) imply that e\ = e 2 = 0. These subcases are already discussed and 
the integrals are equivalent to the integral in equation (|35|) . The only case remaining 
is when 7 = (3 = that gives b\ + b 2 = 0, provided ai 2 a 2 i 7^ 0. We then obtain the 
following first integral: 

H = e (a 21X1 -a 12 x 2 )/b 2 ( ei + ai2XlX2 ^ 

or, equivalently, 

H = a 2 iXi - a u x 2 + b 2 \n(e l + a 12 xix 2 ). (42) 
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2.5. Further notes regarding the known integrals of two-dimensional LV systems and 
quadratic systems 

Many attempts have been made to study the integrability of two-dimensional LV systems 
and general quadratic systems. Different first integrals were found using different 
methods, some similar to those found in this discussion. 

The first integral (126]) along with its conditions f l25|) was probably first found by 
Frommer in 1934 (see Artes and Llibre [21] )• It was also derived by Cairo et ah, who 
used the Hamiltonian method and by Hua et ah, who studied the connection between 
the existence of a first integral and the Painleve property in a general quadratic system. 
In the latter work, the form of the first integral and the vector field are different, but 
through some invertible transformations it is not hard to check that the result is actually 
equivalent. 

Many first integrals were known for the case e\ = and e<i = 0. For instance, 
(13211351) were derived by Nutku [17], Plank [22] and Cairo et al. [3l9l[25] using various 
methods. We remark that the first integral (1351) that has constraints an = 022 = was 
first derived by Volterra himself as a constant of motion (see the book by Hofbauer and 
Sigmund pQ). To our best knowledge, the other first integrals presented here are new. 



3. Three-dimensional LV systems with constant terms 

We consider the following three-dimensional LV systems with constant terms: 
x\ = fi(xi,x 2 ,x 3 ) = £i(&i + a n xi + a 12 x 2 + a 13 x 3 ) + e x , 
x\ = fz(xu %2, x 3 ) = x 2 {b 2 + a 2 ix 1 + a 22 x 2 + a 23 x 3 ) + e 2 , 
x 3 = fa{xi,x 2i x 3 ) = £3(63 + a 3 iXi + a 32 x 2 + 033^3) + e 3 , (43) 

where 6j, (i,j = 1, 2, 3) are arbitrary parameters and (i = 1, 2, 3) are the constant 
terms. In this section, in order to find integrals of the system above we shall make the 
following two Ansatzs for the skew- symmetric matrix T, 

/ -a' -f3' \ 



T 1 (x 1 ,x 2 ,x 3 ) = R 



a' 



\P' 7 



-a 

-i 




(44) 



/ 



respectively, 



T 2 (xi,x 2 ,x 3 ) = R 



( —ax 3 —(3x 2 \ 
—7x1 
7x1 



(45) 



where a, a', (3, /3', 7, 7' G 
function R, namely R = 



ax 3 
\ (3x 2 

are arbitrary parameters. We also use an ansatz for the 



R(xi,x 2 ,x 3 ) = x 



h-lj-z 

Jury 



~ x x l n 1 , where the L 



1, 2, 3) are 



free parameters that are to be determined later on. As discussed in the introduction, 
this Ansatz is restrictive. However, a more general treatise falls outside the scope of the 
present paper. 
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0. 



(46) 



The matrices Ti(xx,x 2 , x 3 ) (i = 1, 2) are integrating factors if and only if curl(T,J) = 
, where /= (fx, f 2 , /a) T - In the three-dimensional case, this condition is equivalent to 

i j k 

d/dxi d/dx 2 d/dx 3 
dHi j dx 1 dH i /dx 2 dHi /dx 3 

where 

VHi = Tif = 1,2). (47) 

In order to find the first integral H\, we expand the above expression with respect to 
the matrix 7\ as follows, 
8H X 



dx\ 
8H X 

dx 2 
dx 3 



-Ra'f 2 - Rf3'f 3 , 
Rot'h - Rj'f 3 , 
R(3'h + Rif 2 . 



(48) 
(49) 
(50) 



and the associated first integral Hx(xx, x 2 , x 3 ) is given by 

Hx(xx,x 2 ,x 3 ) = J (R/3'f 1 + Rj'f 2 )dx 3 + h(x 1 ,x 2 ), (51) 

where h(xx,x 2 ) is found by imposing (j4"8~|) and (1491) . The computation of H 2 is completely 
analogous. 

In the following, we shall derive integrals for the cases ei, e 2 , e 3 7^ 0; ex, e 2 7^ 0, e 3 = 
0; ex 7^ 0, e 2 = e 3 = 0; and ex = e 2 = e 3 = 0. 



3.1. The case e 1; e 2 , e 3 7^ 

A bit of algebra shows that if ex,e 2 ,e 3 7^ 0, the integrating factor matrix T 2 does not 
yield any solutions, so, in this case, we will only use T x . We substitute expressions 
1811501) in the condition (H6l) and obtain 



(52) 



/ (> (Rpfx + Rj'h)- ^-(Ra'fx-Rj'fs) X 



dx 2 
d , 



-Ra'f 2 - Rfi'fs 



dx 3 
d 



dxx 



(R(3'fx + Rif 2 



dx 3 

9 -(Ra'fx ~ Rih) ~ -i-{-Ra'f 2 - Rpfa 



J 



\ dxx dx 2 
Finally, we substitute the vector field (l4"5j) and the function R into the above vector. 
For the first component we find 

-— -\(l 2 - l)6i + (l 2 - l)a n xi + l 2 a 12 x 2 + (l 2 - l)a 13 x 3 + (h ~ 1)— 1 + 

X 2 Xx 

7 [l 2 b 2 + l 2 a 21 xx + (l 2 + l)a 22 x 2 + l 2 a 23 x 3 + (l 2 - 1)— ] - 

x 2 

[(l 3 - l)6i + (l 3 - l)axxXx + (l 3 - l)a 12 x 2 + l 3 a 13 x 3 + (l 3 - 1)—] + 

x 3 Xx 
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7 [^3 + ^3«3i^i + ^32X2 + (h + 1)0332:3 + (h - 1)— ] = (53) 

£ 3 

and the other two components yield similar equations. We now want to find conditions 
on the parameters (a', (3', 7', k, aij, 6j, ti) such that a solution exists. The results for this 
case are summarized in the following lemma. 

Lemma 3.1. The vector field (f73j ) with ei, e 2 , e 3 7^ has a first integral in the following 
cases: 

(i) if the conditions 61 + 62 = 0, 2an + 021 = 0, 2a 2 2 + 012 = 0, a i3 = 023 = are 
satisfied, then the integral is given by 

H = b\X\X2 + a u xlx 2 — a 22 X\x\ + eix 2 — e 2 Xi, (54) 
(ii) if the conditions bi+b 2 = 0, 61+63 = 0, 2an+a 21 = 0, 2an+a 31 = 0, 2a 2 2+«i2 = 0, 
2033 + a i3 = ; and 012013 + 012^23 + ai3«32 = are satisfied then the integral is 
given by 

H = - 261033X1X3 - 2610,22X1X2 - 2aiia 3 3XiX 3 + 1a\\a 22 x\x 2 
+ 2033X1X3 + 2a\ 2 x\x\ + 4022^33X1X2X3 

+ 2(e 3 a 33 + e 2 a 22 )xi - 2eia 33 x 3 + 2eia 22 X2, (55) 

(Hi) if the conditions 6j = and for i 7^ j and i,j = 1,2,3 are satisfied, 

then the integral is given by 

U _ 2 2 2 2 2 2, 2 2,2 2 

n — a lx a 22 x x x 2 — ana 3 3XiX3 — a u a 22 xix 2 -+- aiia 33 XiX3 + a 22 a 3 sX 2 X3 

- ^22033X2X3 + (-aiia 2 2e 2 + aiia 3 3e 3 )xi + (a n a 22 ei - a 2 2a33e 3 )x2 

+ (022^3362 - aiia33ei)x 3 . (56) 

All other first integrals that can be found using Ansatz fiJ4\ ) are related to these three 
cases through a permutation of the coordinates Xj and coefficients {a^, 6j, e^}. 

Proof. First considering terms proportional to xf/xj for i ^ j we find that l\ = l 2 = 
Z 3 = 1. Consequently, we get the following conditions on the parameters, 

«'(6i + 6 2 ) = 0, (57) 

+ 63) = 0, (58) 

7 ' (6 2 + 63) = 0, (59) 

a'(2an + a 21 ) = 0, (60) 

a'(2a 22 + a 12 ) = 0, (61) 

/3'(2an + a 3 i) = 0, (62) 

/3'(2a 33 + a 13 ) = 0, (63) 

7 '(2a 22 + a 32 ) = 0, (64) 

7'(2a 3 3 + a 23 ) = 0, (65) 

-a'ai 3 + /3'ai2 + 7'(a 2 i + a 3 i) = 0, (66) 

tt'^23 + /3'( a 12 + «32) + 7' a 21 = 0, (67) 

a'(ai3 + a 23 ) + f3'a 32 - j'a 31 = 0. (68) 
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(i) We start with the case where a' 7^ 0, 0' = 7' = 0. From (I57H68P the following 
conditions immediately apply: 

61 + b 2 = 0, 2a n + a 2 i = 0, (69) 
2a 22 + ai2 = 0, 013 = ^23 = 0, (70) 
and using floTl) we obtain the first integral (|54p . 

(ii) We turn to the case where a', (3 1 7^ 0, 7' = 0. From (15711651) . the following conditions 
immediately apply, 

61 + b 2 = 0, 2a n + a 2 i = 0, 2an + a 3 i = 0, (71) 

h + 63 = 0, 2a 22 + a i2 = 0, 2a 33 + a 13 = 0. (72) 

While, a' and /3' can be computed from the following linear homogeneous equations 
due to flUSD, 

a'a^ - /3'a 12 = 0, (73) 
a'a 23 + f3'(a 12 + a 32 ) = 0, (74) 
a'(a 13 + a 23 ) + (3'a 32 = 0. (75) 

The solvability condition of the linear system above is given by, 

ai20i3 + oi 2 a 23 + a 13 a 32 = 0. (76) 

If the above condition is satisfied, we obtain the first integral fl55l) due to fl5T|) . 

(iii) Finally we consider the case a', /3', 7' 7^ 0. From the equations (I57H65I) . we have 



61 + ^2 = 0, 2an + a 2 i = 0, 2a 22 + a 12 = 0, 
61 + b 3 = 0, 2a a + a 3 i = 0, 2a 33 + a 13 = 0, 
b 2 + b 3 = 0, 2a 22 + a 32 = 0, 2a 33 + a 23 = 0. 
We simplify the equations above to 



bi = and a,- 



2ajj, 



1 7^ 3i = 1,2,3). 



(77) 
(78) 
(79) 

(80) 
giving 



The parameters a', /3' and 7' can be computed from the equations (1661168 
us the following homogenous linear system, 

— ai 3 a\ 2 (a 2 i + a 3 i) 

«23 (a i2 + a 32 ) a 21 (3 =0. (81) 

y (a 13 + a 23 ) a 32 -a 31 J \ j j 

Substituting (|80|) to the above linear system, we have the following solutions, 

a = aua 22 fi, = -a n a 33 jj, } 7 = a 22 a 33 /i. (82) 

Then if all the parameters of the Lotka-Volterra systems fj43|) satisfy conditions 
above, the system (1431) admits the first integral ( |56l) . 

We have ordered the solutions according to a' 7^ 0, a', (3' 7^ and a', /?', 7' 7^ 0. A 
straightforward computation shows that any other combination yields one of the first 
integrals ( 15411561) with the coordinates and coefficients permuted. For instance, if we 
take a' = 7' = and 0' 7^ 0, we find the first integral (|54l) with the indices permuted 
according to {1,2,3} {1,3,2}. □ 
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We remark that the first integral given in Lemma 13.11 point 1 is the same integral 
that is given in the two-dimensional case. This can be guessed as the first two equations 
of ( 143]) are independent of x 3 . 

In the next three subsections, we now turn to the case where there is at least one 
constant term equal to zero. As it turns out, the first integrals found with Ansatz ( )44l ) 
can all be obtained as limits of the results given in Lemma 13.11 New first integrals can 

83) 
84) 
85) 
86) 
87) 
88) 
89) 
90) 
91) 
92) 
93) 
94) 
95) 
96) 
97) 





00110111011 


R 7 _L R 7 

D\l\ + -fc> 2 i 2 


n 

— u, 


R / 1 R 7 


- u, 


R 7 R 7 


n 

— u, 


A 1 J- A 1 

+ A 23 L 2 


n 


A 7 1 A 7 
/i 32 <i + /i 22 i 3 


= 0, 


^31^2 — ^11^3 


= 0, 


A n h + A 21 l 2 


= -An, 


A 12 h + A 22 l 2 


= ~A 22 , 


A 3 ih + A 21 l 3 


= -^31, 


A 33 l\ + ^23^3 


= -^23, 


A 32 l 2 — A\ 2 l 3 


= ~A 32 , 


A 33 l 2 — A\ 3 l 3 


= A13, 


7e 2 (/ 2 - l)/^2 


= 7e3(^3 - 


/3ei(/i - l)/xi 


= 0e 3 (h ~ 


aei{li — l)/xi 


= ae 2 (l 2 - 



l)/^3 = (^2 - al 3 )e 1 /x 1 = 0, 
l)/x 3 = (al 3 + r yh)e 2 /x 2 = 0, 
l)/x 2 = (/9Z 2 - jh)e 3 /x 3 = 0. 
Here, we have introduced 

Bi = bia - 637 A u = a u oc - a 3i ^ 
B 2 = b 2 a + b 3 (3 A 2i = a 2i a + a 3i (3 

B 3 = h/3 + 6 27 A 3l = a u + a 2i <y. (98) 

Note that 

B l p + B 2l -B 3 a =0 
A U P + A 2i j - A 3i a = Vi. 

3.2. The case ei, e 2 7^ 0, e 3 = 

Given that ei,e 2 7^ 0, e 3 = 0, we conclude that Z x = l 2 — 1 due to conditions f)95H97p . 
and Z 3 must satisfy the following equations: 

- al 3 = al 3 + 7 = 0. (99) 

The results are described in the following lemma. We omit the proof, which is 
straightforward and analogous to that of Lemma 13.11 



An integrating factor matrix method to find first integrals 13 
Lemma 3.2. The vector field (Jty) with ex, e 2 7^ 0, e 3 = has the following first integrals: 



(i) H = biX\X2 + aux\x2 — a 22 X\x\ + e\X 2 — c 2 X\, under conditions 
bi + b 2 = 2a u + a 2 i = 2a 22 + a 12 = and a 13 = a 23 = 0; 

(ii) H = (013 — a 23 )x\x 2 x\j2 + e\x 2 x 3 — e 2 x\x 3 , under conditions 
Oi3 - ^23 ^0, b 1 + b 3 = b 2 + b 3 = 0, a 2 i + a 31 = a u - a 21 = 0, 
«22 + 032 = «i2 - «32 = 0, and a 13 + a 23 + 2a 33 = 0; 

(Hi) H = [e\X 2 — e 2 Xi)x 3 3 , where l 3 is a solution of 

h - b 3 l 3 = a n + a 31 l 3 = a 12 + a 32 l 3 = a 13 + a 33 l 3 = 0. 

under conditions b\ + a 31 + a§ 2 + a§ 3 ^ 0, b\ — b 2 = 0, an — a 2 i = 
and b 3 au — b\a 3 i = for i = 1,2,3 

These are all first integrals that can be found using Ansatz J73p . 

We also remark that the first integral given in Lemma 13.21 point 1 is trivial as it 
follows directly from the two-dimensional case. 

3.3. The case e\ ^ 0, e 2 = e 3 = 

In this case, it immediately follows that h = 1 due to conditions (|96|) and (1971) . and we 
have also the following equation to satisfy condition ( J95l) . 

(3l 2 - al 3 = 0. (100) 

All results for this case are given in the following lemma, 



Lemma 3.3. The vector field with e\ ^ 0, e 2 = e 3 = has the following first 
integrals: 

1. H = —B 2 X\ — A 2 ix\/2 + aei In \x 2 \ + (5e\ In \x 3 \, where a, (3 are solutions of A 22 = 0, 
A 23 = 0, under conditions 

bi = an = a 12 = a 13 = and a 22 a 33 - a 23 a 32 = 0; 

2. H = —a 23 X\X 3 + eiln \x 2 \, under conditions 

&2 = «2i = 0-22 = 0, a 23 ^ and b x + b 3 = a u + a 3i = for i = 1,2, 3; 

3. H = pin \x 3 \ + a In \x 2 \, where a, (3 are solutions of B 2 = and A 2i = 0, 

for i = 1, 2, 3 under the condition (b 2 , a 2 i, a 22 , a 23 ) T = X(b 3 , a 3 ±, a 32 , a 33 ) T for some 
A G R; 

4- H — A 33 x\x 3 + (3ei In \x 3 \ — 7eiln|x 2 | ; where 7 are solutions of B 3 = 0, A 3 i = 0, 
A 32 = 0, under conditions 

A 33 ^ 0, h + 6 3 = an + a 3i = for i = 1, 2, 3 and (61, a n , an) T = A(6 2 , a 21 , a 22 ) for 
some A G M; 

5. H = (a 13 + a 23 )a;3 — (a 12 + a 32 )x 2 + e\ In \x 3 \ — e\ In |x 2 |, under conditions 

ai2 + a 32 7^ 0, a i3 + a 23 ^ 0, bi + 6 2 = bi + 63 = 0, an + a 2i = an + a 31 = and 
«i2 + 022 = ai 3 + a 33 = 0; 
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6. H = —023X1X2 £3 + eix 2 2 /I2, where I2 = — (013 + a 33 )/a 2 3, under conditions 

a 23 7^ 0, a V 3 + a 33 7^ 0, b x = a 2i = a 22 = and b x + o 3 = a u + a 31 = a 12 + a 32 = 0; 

7. H = (bi + anXi)xiX 2 2 x 3 3 + eix 2 2 x 3 3 , where l 2 ,l s are solutions of b 2 h + W3 — — °i; 
o 2 i/ 2 + 031Z3 = — 2on ; o 22 Z 2 + o 32 / 3 = o 23 Z 2 + o 33 / 3 = 0, wider conditions 

ai2 = ai3 = 0, a 22 (-6ia3i + 26 3 an) + a 32 (-2& 2 oii + &1O21) = and a 23 (-&ia 3 i + 
263011) + a 33 (-26 2 an + 6 x a 2 i) = 0; 

8. H = (013 + a 3 3)xix 2 2 x 3 3+1 + eia; 2 2 x l 3 , where / 2 , Z 3 are given below 

, _ 7(^13 + 033) , /3(oi3 + o 33 ) 
2 ~ 4 3 ~ A ' 

^-23 ^23 

and /3,7 are solutions of equations B 3 = 0, A31 = 0, A 32 = under the conditions 
A33 ^ 0, A23 7^ 0, a i3 + a 33 ^ 0, 61 + 63 = an + a 3i = a 12 + a 32 = and 
(61, an, a i2 ) T = A(6 2 , O21, o 22 ) T /or some A G M; 

9. H = ((012 + a 22 )x 2 /7 3 + (a i3 + a 23 ):r 3 /(Z 3 + l))xix 2 x 3 3 + eix 2 2 x 3 3 /l 3 , u>/iere / 2 = 
— (012 + o 22 )/(a 22 — a 32 ) ana 1 Z 3 = — £ 2; under conditions 

012 + o 22 7^ 0, a 22 - a 32 7^ ; a i3 + a 33 7^ ; a 23 - a 33 7^ 0, h + 63 = b 2 - 63 = 7 
on + o 3 i = a 2i - a 3 i = and (a 13 + a 33 )(a 22 - a 32 ) - (a i2 + a 22 )(a 23 - a 33 ) = 0. 

We note that the integrals in the lemma 13.31 point 3 and 5 do not depend on the 
variable x\. They immediately follow from the two-dimensional case. 



3. If.. The case e\ = e 2 = e 3 = 

Finally, we shall discuss the case where all the constant terms are zero. The equations 
(1951 - [97|) are satisfied immediately. This means we only need to find conditions on the 
parameters in order to satisfy equations (l83l- l94"l) . 

In the following lemma, we describe our results obtained using an integrating factor 
matrix of the form T 2 . 



Lemma 3.4. The vector field (43) with e\ = e 2 = e3 = has the following first integrals: 

1. H — a(b± In |x 2 | — 6 2 In \xi\ — a 2i xi) + j3(b\ In |x 3 | — 63 In \x±\ — 031X1) where a, (3 are 
solutions of equations A 2 2 = and A23 = 0, under conditions 

61 7^ 0, an = 012 = a i3 = and a 22 a 33 — a 32 a 23 = 0; 

2. H = a(b 2 /xi — a 2 i In \x%\ + an In |x 2 |) + /3(6 3 /xi — a 3 i In \x± \ + an In \x 3 \) where a, (3 
are solutions of equations A 2 2 = 0, A23 = 0, under conditions 

an 7^ 0, 61 = 012 = a i3 = and a 22 a 33 — a 32 a 23 = 0; 

3. H = y4 31 ln|x 3 | + An In |x 2 | — v4 21 ln|xi| + A 22 a; 2 /:ri where a, /?, 7 are solutions of 
B\ = 0, A13 = and /3 + 7 = ; under conditions 

A n + A li 0, ^22 7^ 0, 61 - 6 2 = a 12 - a 22 = a i3 - a 23 = and b x a 33 - 6 3 ai 3 = 0; 

4. H = -(an - a 2i ) In \x 3 \ + (a u - a 3l ) In \x 2 \ - (a 2i - a 31 ) In |xi| + (a 22 - a 32 )x 2 /xi - 
(ai 3 — a 23 )x 3 /xi ; under conditions 

(an - o 3 i) 2 + (an - a 2i ) 2 7^ 0, a 22 - a 32 7^ 0, a 23 - a 33 7^ and b x - 6 3 = b x - 6 2 = 
a 12 — a 22 = 013 — 033 = 0; 
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5. H = xfx 2 , where /3,7 are solutions of B 3 = 0, A 31 = 0, A 32 = and A 33 = 0, 
under conditions 

(bi,au,a 12 , a 13 ) T = X(b 2 , a 21 , a 22 , a 23 ) T for some A G R; 

6. H = x^x 2 +a x 3 , where a,/?, 7 are solutions of B 2 = 0, B 3 = 0, A 2i = and A 3i = 
(% = 1, 2, 3 ), under conditions 

(bi, an, a 12 , a 13 ) T = Xi(b 2 , a 21 , a 22 , a 23 f and 

(b 2 , a 21 , a 22 , a 23 ) T = A 2 (o 3 , a 31 , a 32 , a 33 ) T for some Ai, A 2 G R; 

la. H = x l lx 2 2 (A 12 x 2 /(l 2 + 1) + A 33 x 3 ) where h = -f3(-a 23 + a 33 )/A 33 , l 2 = -A 13 /A 33 
and a, (3, 7 are solutions of B 3 = 0, A 31 = 0, A 32 = and a + [3 = 0, under 
conditions 

A 33 ^ 0, A 13 ^ 0, A 23 ^ 0, A 12 ^ 0, b 2 - b 3 = a 21 - a 31 = and {b 1 ,a 11 ,a 12 ) T = 
X(b 2 , a 2 i, a 22 ) T for some A G R; 

7b. H = x l 1 x l 2 2 (-b 3 (3/l 1 + A 33 x 3 ) where l x = -a 33 (3/A 33 , l 2 = -(a 33 j)/A 33 and (3, 7 are 

solutions of B 3 = 0, A 3i = 0, A 32 = 0, under conditions 

h 7^ 0, a 33 7^ 0, A 33 7^ 0, a 31 = a 32 = 0, and 

(61, an, ai 2 ) T = A(o 2 , a 2 i, a 22 ) T for some A G R; 
7c. H = x' 1 1 x 2 2 ((a 2 i - a 3 i)x 1 /(l 1 + 1) + (a 22 - a 32 )x 2 /l 2 + (a 13 - a 23 )x 3 ) where h = 

(a 23 - a 33 )/ (ai 3 - a 23 ) 7 / 2 = (a 33 - ai 3 )/ (a i3 - a 23 ) ; under conditions 

au ~ a 3 i 7^ 0, a i2 - a 32 7^ 0, a i3 - a 33 7^ 0, a 23 - a 33 7^ 0, a 13 - a 23 7^ 0, 

bi — b 2 = b\ — b 3 = and an — a 2 i = ai 2 — a 22 = 0; 

7c?. H = x l ^x l 2 (-(b 2 + b 3 )/(li) - (a 2i + a 31 )x x /l 2 + (a 13 + a 23 )a; 3 ) 7 w/iere k = -(a 23 + 
^33)/ ( a i3 + ^23) «™d ^2 = -(033 - a 13 )(a 13 + a 23 ), under conditions 
ai3 + «23 7^ 0, 61 + b 2 = 0, an + a 21 = 0, a 12 + a 22 = 0, a 12 - a 32 = and 
(61 - 6 3 )a 23 - (6 2 + 03)^13 = 0; 

8a. H = x l fx%x l i(bi + onii) where k = -(a u B 2 ) / '(a u B 2 - M21), l 2 = -(b 1 a)/(B 2 ), 
h — (bi(3)/(B 2 ), and a, (3 are solutions of A 22 = 0, A 23 = 0, under conditions 
b\ + a n 7^ 0, B 2 7^ 0, an-B 2 - b x A 2X 7^ ; a 12 = a 13 = and a 22 a 33 - a 32 a 23 = 0; 

8b. H = x l ^x l 2 x l 3 (ai 2 x 2 /l 3 — a\ 3 x 3 jl 2 ) where 

l _ (a 22 - a 32 )(a 23 - a 33 ) ^ 012(^23 ~ ^33) 

-ai 2 (a 23 - a 33 ) + a 13 (a 22 - a 32 ) -ai 2 (a 23 - a 33 ) + ai 3 (a 22 - a 32 ) ' 

and l 2 + l 3 = —1, under conditions 

ai3 7^ 0, a 23 7^ 0, a 22 -a 32 7^ ; a 23 -a 33 7^ 0, 61 = an = andb 2 -b 3 = a 21 -a 31 = 0; 

8c. H = x l iX l 2X l £((bi + b 2 )/(l 3 ) + (an + a 2 i)x 1 /l 3 + (a i2 + a 22 )x 2 /l 3 - (a 23 + a 33 )x 3 /l 1 ) 
where 

_ A 22 (A 21 - A n ) _ A n (A 12 -A 22 ) A 31 (A 33 - A 23 ) 

A\\A 22 — A 2 \A\ 2 A\\A 22 — A 2 \A\ 2 A 31 A 23 — A 2 \A 33 

where A\\A 22 — A 2 \A\ 2 7^ ; A 3 \A 23 — A 2 \A 33 7^ and a = (3 = 7 under conditions 
B 1 A 2 2(A2i-A 11 )+B 2 Au(A 1 2-A22) = 0, A 13 A 22 (A 21 - A n ) , +A 23 A U (A 12 -A 22 ) = 
0, B 3 A 23 (A 21 -A 31 ) + B 2 A 31 (A 33 -A 23 ) = and A 32 A 23 (A 21 - A 31 ) + A 22 A 31 (A 33 - 
A 23 ) = with a — j3 — 7. 
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Remark that the integral in the Lemma [3.41 no 5. is the same integral we found in 
the two-dimensional LV case. 

3.5. On the first integrals of three-dimensional LV systems 

We here give some remarks about lemma 13. 4[ in which we discussed first integrals 
of three-dimensional LV systems that have already been extensively discussed in the 
literature [71[Tnl[TT|[T3HT6||2T||22|[25|[26] . The forms of first integrals obtained in lemma 
13.41 (points 1-4) are similar to the ones obtained by Plank (1995) |22J . which in fact, are 
special cases [27] of invariants found in Cairo and Feix (25]. It is not hard to show that 
their integrals are different from the ones we have obtained in this paper. 

The other integrals in the lemma l3~41 (points 5-8) have the following general form 

xl'x^xl 3 (f{x 1 ,x 2 ,x 3 ) Xi (101) 

where Ai, A 2 , A 3 and A 4 are some constants and if is a polynomial function of degree 
one in xi,X2,xa. To our best knowledge, the above general form first appeared as 
a first integral in Cairo and Feix [25], except the fact that the first integral in their 
paper is time-dependent. Through a time-rescaling (see Hua et al. [20J) they obtained 
a time-independent first integral. The existence of first integrals of this form of three- 
dimensional LV systems is also extensively investigated by Cairo [ID] . He investigated 
a polynomial function tp of degree one and two. The integral functions in point 5-6 has 
the same form with the ones he found [IDl Theorem 2(1)]. The forms of integrals in 
point 7 do not seem to have been recognized before, thus these new results extend the 
known results on integrals of the form (jlOip . The integral of the form 8a generalizes 
integrals obtained in [TQl Theorem 2(8-13)]. Finally, integrals of the form 8b and 8c 
seem to be new. 

4. Conclusion 

In this paper, we have derived first integrals for two and three-dimensional LV systems 
with constant terms through an integrating factor matrix. We make Ansatzs and obtain 
conditions for the existence of first integrals. By our method, the search for integrals in 
dynamical systems changes into a linear algebra problem. We note that some conditions 
of the two and three LV system with constant terms to have first integrals do not involve 
any constant terms. This is good because under such conditions tfirst integrals are 
preserved for any constant terms. 

In the two-dimensional case, the integrating factor matrix T(x\,X2) in equation fl4]), 
along with the condition such that T is an integrating factor, turns out to be similar 
to the one, used by Plank [22] due to S — T -1 . In the three-dimensional case this 
property no longer applies, clS 3j 3 X 3 skew-symmetric matrix is not invertible. Hence, 
our analysis generalises the work of Plank's. 

Comparison with Darboux method The Darboux method has been applied to 
find first integrals of two-dimensional [HIE] and three-dimensional [TUIEH] LV systems. 
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Compared to this method, our method has advantages in the context of searching for 
a first integral of LV systems with constant terms. To apply the Darboux method, one 
must seek an algebraic curve of a vector field, and LV systems without constant terms 
have natural algebraic curves as the axes are invariant. For systems with constant terms, 
this no longer applies. 

Comparison with Hamiltonian method The existence of first integrals of two- 
dimensional systems has been obtained using Hamiltonian method [T8H20] . They 
assumed that the integrals are products of two or three polynomial functions of degree 
one. The advantage of our method is that we only make a single Ansatz. Gao [15] has 
also applied a Hamiltonian method to find first and second integrals of special cases of 
three-dimensional LV systems, where the linear terms are absent. 

Comparison with Frobenius method The Frobenius method was first introduced by 
Strelcyn and Wojciechowski [12] to find a first integral for three-dimensional systems. It 
has been used to find integrals for LV systems by Grammaticos et al. [11] . Unfortunately, 
they only look at a special case of the LV system, which is the so-called ABC system. 
Comparison with Carleman embedding method The existence of first integrals 
in n-dimensions has been studied [HI [25] through the Carleman embedding method. 
However, the integrals that are obtained are time-dependent and are refered to as 
invariants. 

For future work, it would be a challenging problem to find a more general integrating 
factor matrix of a vector field in dimension greater than or equal to two. Several first 
integrals with the corresponding conditions have been published in the literature and 
we could use these results to endeavour to find the general integrating factor matrix. 
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